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Abstract 

A point P 6 M™ is represented in Parallel Coordinates by a polyg- 
onal line P (see |Ins99| for a recent survey). Earlier Ins85 j, a surface 
a was represented as the envelope of the polygonal lines representing 
it's points. This is ambiguous in the sense that different surfaces can 
provide the same envelopes. Here the ambiguity is eliminated by con- 
sidering the surface a as the envelope of it's tangent planes and in 
turn, representing each of these planes by n-1 points |Ins99| . This, 
with some future extension, can yield a new and unambiguous rep- 
resentation, a, of the surface consisting of n-1 planar regions whose 
properties correspond lead to the recognition of the surfaces' proper- 
ties i.e. developable, ruled etc. Hun92 ) and classification criteria. 

It is further shown that the image (i.e. representation) of an al- 
gebraic surface of degree 2 in M." is a region whose boundary is also 
an algebraic curve of degree 2. This includes some non- convex sur- 
faces which with the previous ambiguous representation could not be 
treated. An efficient construction algorithm for the representation of 
the quadratic surfaces (given either by explicit or implicit equation) 
is provided. The results obtained are suitable for applications, to be 
presented in a future paper, and in particular for the approximation of 
complex surfaces based on their planar images. An additional benefit 
is the elimination of the "over-plotting" problem i.e. the "bunching" 
of polygonal lines which often obscure part of the parallel-coordinate 
display. 
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1 INTRODUCTION 

Our purpose here is expository, sparing the reader from most of the mathe- 
matical tribulations and, focusing on the more intuitive aspects of the repre- 
sentational results. After a short review of the fundamentals, the essentials 
of the mathematical development are given together with some detailed ex- 
amples to clarify the nuances and satisfy the more mathematically inclined. 

In parallel coordinates (abbr. ||-coords), a point in R 2 is represented by a 
line and a line is represented by a point yielding a fundamental point <-> line 
duality. There follows the representation of p-flats (planes of dimension 

2 < p < n — l)in M. n in terms of indexed points jlns85j . Naturally, for 
non-linear objects the representation is more complex, especially if they are 
also non-convex. The points of a curve in R 2 can be mapped directly into a 
family of lines whose envelope defines a curve ( "line-curve" ) . Actually this 
is awkward and also clutters the display. Instead we map the tangents of 
the original curve into points to obtain the "point-curve", sometimes called 
"dual-curve", image directly as shown in Fig. ^ In short, this approach 
provides a convenient point-to-point mapping Ins99 . 




Figure 1: Point-curve mapped into point-curve, result of considering the 
point-curve as the envelope of it's tangents. 

Applying these considerations it was proved that the image of an alge- 
braic curve of degree n is also algebraic of degree at most n(n — 1) in the 
absence of singular points jlzhOlj . This theorem is a generalization of the 
known result that conies are mapped into conies Dim84 in six different 
ways. 




Figure 2: For a hyper-surface, it's representation's boundary curves are 
determinedfrorn the images of curves contained in hyper-surface. 



Perhaps we are "pushing our luck" , our intent is to apply next the point- 
to-point mapping in the representation of surfaces considered as the envelope 
of their tangent planes; with the resulting image being constructed from the 
representation of tangent planes |Ins99j . As has already been pointed out, 
planes can be represented in ||-coords by indexed points. The collection 
of these planar points, grouped for each index, is the representation of the 
surface. 

In the past |Ins85j . surfaces were represented in ||-coords by the envelope 
of the polygonal lines representing the surfaces' points. By itself this is 
ambiguous. For example the image of a sphere in n-dimensions is the same as 
the image of the surface obtained by the intersection of n cylinders properly 
aligned having the same radius. In applications this was ameliorated by 
accessing the correct equation of the surface when needed. Not only is 
the ambiguity completely removed with the new representation, but also 
non-convex surfaces can be nicely treated something that was not possible 
previously. 

Hung Hun92 first applied this notion and found that regions represent- 
ing developable surfaces consists only of the boundary curves (i.e. there are 
no interior points), and also that ruled surfaces can be recognized from char- 
acteristic properties of their corresponding regions. Encouraged by these 
initial results the analysis is extended to more general surfaces yielding use- 
ful criteria in the approximation of surfaces by simpler ones; but we are 



getting ahead of ourselves. 

At first we lay the foundations, then derive the representation of quadratic 
algebraic surfaces and further generalize to higher dimension as well as more 
complex hyper-surfaces. As a result an efficient algorithm for constructing 
the representation of the quadratic surfaces (given either by explicit or im- 
plicit equation), and a proof that the image of an algebraic surface of degree 
2 in W 1 is also an algebraic curve and of degree 2 are obtained. 

2 GENERAL REPRESENTATION OF HYPER- 
SURFACE 

Tn general, the method employed below applies to the class of smooth hyper- 
surfaces in M. n having a unique tangent hyper-plane at each point. Equiva- 
lently, each such hyper-surface is the envelope of it's tangent hyper-planes. 
This is our point of departure, for it enables us to represent each tangent 
hyper-plane in ||-coords by n — 1 indexed points Ins!)!) . The hyper-surface's 
representation consists of the n — 1 points sets, one for each index [Hun92 . 
For the present we restrict our attention to algebraic hyper-surfaces and 
in particular those defined by quadratic polynomials. To simplify matters, 
most of the analysis is done in 3-dimensional space but in a way which points 
to the generalization for M. n . 

2.1 Hyper-Planes Representation 

An n-dimensional hyper-plane tt in W 1 

n 

tt : y^CjXj = cq, (1) 
1=1 

is represented by the n — 1 indexed points Ins99 . For our purposes only 
the first 

n n 

7fi...n = ~ 1)q,c ,^q). (2) 

1=1 i=l 

needs to be studied. The remaining n — 2 points have similar form differing 
only in the factor (i — 1) of the c%. An important property is that the 
horizontal distance between the i-adjacent (in the indexing) points is the 
equal to the coefficient q; from which the sequence of indexed points can be 
generated from the coefficients or vice-versa. 



2.2 Hyper-Surfaces Representation 



Moving on to the representation of non-linear hyper-surfaces in WL n from 
their tangent hyper-planes. Let a be a smooth ra-dimensional hyper-surface 
generated by the differentiable function F(xi, . . . , x n ) = 0, and an arbitrary 
point (x\, . . . , x„) G a. Then the hyper-surface's tangent hyper-plane at this 
point is given by : 

n Pit? 

J>-z°)^(z?,...,4) = 0. 
i=i 

Taking the x® as parameters, the coefficients of the tangent hyper-planes 
can written as a function of points which satisfy the hyper-surface's equa- 
tion. Namely, the family of tangent hyper-planes of a is represented in 
homogenous coordinates, by a collection of sets a (see Fig. |2j), containing 
the indexed points representing each member (i.e. hyper-plane) of the fam- 
ily. Each of the n — 1 indexed set, o~j... n i'...(j—iy^ consists of the points with 
the same index. 

In the remainder the analysis is confined to the first indexed set ai... n 
using a shorter notation a defined as, 

* = {(P(x), S(x), Q(x)) \xea)}, (3) 

where for o\,„ n and the tuple x = [x\, . . . , x n ) : 



O 771 

i=i 1 

«*> - ± 8 £,- 

i = l 

In general, the representation in || -coords is constructed via the rational 
transformations 

_ P(x) _ 5(a) 

The representation of all these hyper-planes transform an n-dimensional 
hyper-surface into subsets of P 2 ; regions which are distinguished from each 
other by their indices. The algorithm which constructs and describes these 
regions is presented next. 




Figure 3: The boundary of each region is an image of a curve embedded in 
surface. 

3 REPRESENTATION OF QUADRAT- IC HYPER- 
SURFACE 

At first we treat the class of algebraic surfaces in M 3 described by quadratic 
polynomials. Mercifully, the corresponding system of transformations P|) 
can be linearized. The next step is to determine the boundary of the regions 
representing the surface. Without getting into details, the existence of the 
boundary can be assured by selecting an appropriate spacing of axes in the 
system of the || -coords, which eventually reflect by changing the constant 
multipliers of the first equation in J3J) . 

3.1 Definition of the Regions' Boundary 

Let a be a quadratic surface whose representation is the region a. The 
boundary points are those whose every neighborhood contains both inte- 
rior and exterior points. For this case both the transformation @ and the 
surface, are defined by polynomials and hence are differentiable. The basic 
properties including continuity are therefore preserved under the transfor- 
mation. 

Geometrically we rely on the differentiability in finding those points p € 
a so that we can "move" from p in any direction and still remain in region; 
these are interior points of a. Clearly the boundary points are easily found as 
the complement of the interior of a. The condition for determining whether 



a point is interior or not is given by theorem of implicit function Mar85 . 
Equivalently, a point a £ a is interior point if and only if the Jacobian, 
J(F)\a, at this point is different from zero. Conversely, a point b € a for 
which this is not true is necessarily a boundary point; namely, a point b E a 
s.t. J(F)\i = 0. In essence the theorem tells us that a is closed set and the 
complement of its interior is the sought after boundary. 

Generalizing the above for R n we get the mapping [x\ } . . . , x n ) — > (77, £, -0) 
into the projective space, were x = ^ and y = 4. Restating the condition 
in terms of differential products using homogeneous coordinates with the 
variables 77, £ and ip yields, 

(•qd^dip - £dqd\l) + ^dr]d()dF = 0. (6) 

This form is more convenient for handling hyper-surfaces embedded in W n 
where n > 3, for n = 3 eq. (jHJ) can be written equivalently in term of 
Jacobian as, 

J(F)ip 3 = Det 
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Substituting 77, £ and t/> in terms of the Xj's yields an equation which defines 
an algebraic surface, a', in M 3 . Geometrically, the boundary 7 consist of 
points which represent tangent planes touching at points of a on the inter- 
section a 1 fl c Hence, 7 is the image of the algebraic curve a 1 H a = 7 (see 
Fig. EJ). 

Combining the criterion, eq. 0, for the boundary with the equation 
of the surface (embedded in 3-dimensional space) and the transformation 
equations (in homogenous coordinates) yields : 

F(x 1 ,X 2 ,X 3 ) = 0, 

r]Q(x 1 ,x 2 ,x 3 ) - TpP(xi,x 2 ,x 3 ) = 0, , > 

£Q(x 1 ,x 2 ,x 3 ) - ipS(xi,x 2 ,x 3 ) =0, 

J(F(xi,x 2 ,x 3 ))V 3 =0. 

Solving for 77, £ and ip yields the equation of the boundary. Note that if F 
is a polynomial of degree 2, then the degree of J(F)xjj 3 is < 2 in terms of all 
variables, while it is linear in terms of x\, x 2 and x 3 . 

Thus far we have constructed a system of four equations (jSJ in six vari- 
ables which define a mapping from the M 3 into the projective plane P 2 . Our 





aim, is to determine the specific equation of the region's boundary explic- 
itly. This involves solving this system of equations in terms of rj, £ and i/j 
by eliminating the variables x\, X2 and X3. 

The equation's structure turns out to be very advantageous. Since the 
last three equations are linear, the elimination can be done by isolating 
a variable (finding an explicit expression in term of the other variables), 
and substituting in the remaining linear equations. When all is said and 
done, each of the variables xi, X2 and X3 can be expressed as a rational 
equation in rj, £ and ip. Upon substitution of these expressions into F the 
boundary's equation in homogeneous coordinates is obtained. It follows that 
the boundary 7 is a quadratic curve. 

3.2 Algorithm 

The algorithm's input is an equation of algebraic surface a : F(x±, X2, X2) = 
of degree two and the output is the polynomial which describes the bound- 
ary of the surface's image in ||-coords. It is noteworthy that the algorithm 
applies to implicit or explicit polynomials with or without singular points. 
The formal description is followed by examples which clarify the various 
stages and their nuances. 

For a given polynomial equation F{x\, X2, x%) = of degree 2 and a 
spacing of axes Si... ra : 



• Let : 

• Write the three linear equations: 

«) ^Li(i-i)f-^Li| = o, 

c) J(F)iP 3 = 0. 

• Using substitution write 

Xi = fi{v,£,ip) , for i = 1,2,3. 

• Substitute 

F(fi(v, £, V0> f2(v, 6 1>),f 3 (v, & VO) = o. 

• Retain the equation's numerator. 

• The output is obtained by substitution: 

v <- x > £ <- y » ^ ■<- 1. 

All this falls into place with the following examples. 

4 EXAMPLE OF QUADRATIC SURFACE AND 
THEIR TRANSFORMS 

In the first example is quite detailed to accommodate the readers wishing 
to follow the application of the algorithm in depth. Let a be 3-dimensional 
saddle (see Fig^J generated by the polynomial equation, 

F(x) = (^) 2 - (^) 2 + x 3 = x\ - x\ + 4x 3 = 0, where x = (xi,x 2 ,x 3 ), 
and Si___ n the standard spacing of axes. 

step 1 Let : 

m= £-=i(*-i)lg = -2^2+8, 

S(x)= ^ =iX .g_ 2F = _ AX3) 



Figure 7: Hyperboloid of two sheets a: x 2 — Ay 2 + 2z 2 = —2 — > 7 
x 2 - 2xy + Ay 2 - 1 = 0. 



step 2 Write 1 the three linear equations: 

a) r]Q(x) — ipP(x) = 

2(77x1 + (V> — ??)X2 + 277 — Aip) = 0, 
6) £Q(x) - ipS(x) = 

2(£xi - £x 2 + 2^x3 + 2C) = 0, 
c) J(F(x))tjj 3 = (xi(t/> - 77) + X 2 7? - 0. 

Notice: substitution of 77, £ and tp in terms of x\, x 2 and X3 yields a 
surface in M 3 , 

a' : 32(-x| + 4x 2 + x\ + 2x 3 - 2xj) = 0. 

Hence 7 = a' fl cr, 

cf - x 2 + 4x 3 = 

-x| + 4x 2 + + 2x 3 - 2xi = 

step 3 Using simple substitution write 

b) => ^3 = - ^2; 2+2) , 



lr The authors acknowledge and are grateful for the use of the symbolic manipulation 
program Singular developed by the Algebraic Geometry Group, Department of Mathe- 
matics, University of Kaiserslautern, Germnay. 



Then using equation a) we get, 

_ 2^-^-6^+4^ 
X2 V>(V>-2»?) ' 

_ g(-g+2q+2V>) 
X3 2</)(V>-2»7) • 

step 4 Substitute in i 7 , 

x l ~~ x 2 + ^3 = 

16^ 2 - 16»/)??-4^g+g 2 -iriZ+ir? _ n 

step 5 Retain the equation's numerator, 

16V> 2 - 16V?r? - 4V>£ + - 4??£ + 4r/ 2 = 0. 

step 6 Finally, the output is obtained by substitution, 
rj <- x , i^y , ^ <- 1: 
7 : 16 — 16x — 4y + y 2 — Axy + 4x 2 = 0. 

The surface and its image including the boundary curve are shown in FigHJ 
The representation of other quadratic surfaces is illustrated in figures EJ El 
andEI 



5 CONCLUSION 

The new representation 

• is constructive, 

• enables the representation of non-convex objects, 

• maps algebraic surfaces to regions having algebraic curves as bound- 
aries, 

An important "fringe benefit" is the avoidance of the "over-plotting" prob- 
lem in 1 1 -coords where polygonal lines obscure portions of the display. 
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